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Abstract 

We investigate integrable second order equations of the form 

P{'^xxi '^xyi ^yyi '^xti Uyt, ^tf) — 0. 

Familiar examples include the Boyer-Finley equation Uxx + Uyy = e"" , the potential 
form of the dispersionless Kadomtsev-Petviashvili (dKP) equation Uxt — \u^x — '^yy^ 
the dispersionless Hirota equation (a — f3)e^''y + {P — 7)e"«* + (7 — a)e"*=" = 0, etc. 
The integrability is understood as the existence of infinitely many hydrodynamic 
reductions. We demonstrate that the natural equivalence group of the problem 
is isomorphic to Sp{6), revealing a remarkable correspondence between differential 
equations of the above type and hypersurfaces of the Lagrangian Grassmannian. 
We prove that the moduli space of integrable equations of the dispersionless Hirota 
type is 21-dimensional, and the action of the equivalence group Sp{6) on the moduli 
space has an open orbit. 

MSG: 35Q58, 37K25, 53A40, 53B25, 53Z05. 

Keywords: Integrable equations, Hydrodynamic Reductions, Lagrangian Grass- 
mannian, Hypersurfaces, Rational Normal Gurve, GL(2)-structure. 
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1 Introduction 



We investigate a general class of three-dimensional second order equations of the form 

F{UxX, Uxy, Uyy, Uxt, Uyt, Uu) = (l) 

where u = u{x,y,t) is a function of three independent variables. Equations of this type 
arise in a wide range of applications including non-linear physics, general relativity, differ- 
ential geometry, integrable systems and complex analysis. For instance, the dKP equation, 

also known as the Khokhlov-Zabolotskaya equation, arises in non-linear acoustics [52], as 
well as in the theory of Einstein- Weyl structures [18]. The Boyer-Finley equation, 

^xx ~t~ "^yy ^ 5 

which is descriptive of a class of self-dual 4-manifolds, has been extensively discussed in 
the context of general relativity [12]. The equations 

Hess u = 1 and Hess u = Au, 

where Hess is the determinant of the Hessian matrix of u, and A is the Laplacian, 
appear in differential geometry in the theory of affine spheres and special Lagrangian 
3-folds, respectively [13, 32]. A subclass of equations of the form (1), 

was discussed recently in [41] in connection with hydrodynamic chains satisfying the 
so-called Egorov property. Equations of the form (1) typically arise as the Hirota-type 
relations for various (2-|-l)-dimensional dispersionless hierarchies. For instance, equations 
of the dispersionless Toda hierarchy [49, 14, 53], 

(A - /i)e-^^^^" = Ae-^=-°^" - /ie"^-^'^", 
(A - /i)e^^-^^" = Ae^^-^^" - /xe^^-^^", 

take the form (1) after one replaces the 'vertex' operators Dx, D^, Dx, by ordinary 
partial derivatives. Equations of the dispersionless Toda hierarchy play an important role 
in complex analysis underlying the integrable structure of the Dirichlet problem for simply- 
connected domains [51, 39]. A generalization of this construction to multiply-connected 
domains [34] leads to a remarkable hierarchy of Hirota-like equations which constitute the 
so-called universal Whitham hierarchy [33]. Further examples of Hirota-type relations 
arise in the theory of the associativity (WDVV) equations [8, 10, 30], e.g., 

pUxx+Uxy+Uxt — Uyt pUxy+Uyy+Uyt—Uxt _|_ pUxt+Uyt+Utt—Uxy _ n 
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^utt-uxx _ sinhtij^y/ smhuyti 

coth u.j.y = coth Ur^t coth Uyt, 

etc. These equations arise as certain differential constraints which should be considered 
along with a full set of the WDVV equations. 

Equations of the above type have been approached by a whole variety of modern tech- 
niques including symmetry analysis, differential-geometric and algebro-geometric meth- 
ods, dispersionless 5-dressing, factorization techniques, Virasoro constraints, hydrody- 
namic reductions, etc. However, until recently there was no 'intrinsic' approach which 
would explain the integrability of these and other examples. Moreover, there was no sat- 
isfactory definition of the integrability which would (a) be algorithmically verifiable, (b) 
allow classification results and (c) provide a scheme for the construction of exact solutions. 
We emphasize that equations of the form (1) are not amenable to the inverse scattering 
transform, and require an alternative approach. Such approach, based on the method 
of hydrodynamic reductions and, primarily, the work [26], see also [14, 22, 23, 38], etc, 
was proposed in [20]. It was suggested to define the integrability of a multi-dimensional 
dispersionless system by requiring the existence of 'sufficiently many' hydrodynamic re- 
ductions which provide multi-phase solutions known as non-linear interactions of planar 
simple waves. Technically, one 'decouples' a three-dimensional PDE (1) into a pair of 
commuting n-component (1 + l)-dimensional systems of hydrodynamic type, 

Rl = X'{R) -R^, Rl = Rx, (2) 

where the characteristic speeds A* and satisfy the commutativity conditions 

— A* jj,^ — jj} 

dj = djij , see [50] . The best way to illustrate the method of hydrodynamic reductions is 
to discuss an example. 

Example. Let us consider the dKP equation, Uxt — \u'ix = '^yy^ ^'^d introduce 
the notation Uxx = cl, Uxy = b, Uxt = P, Uyy = p — \a?] this results in the equivalent 
quasilinear representation of the dKP equation, 

ciy = bx, at = px, bt = Py, by = {p- ^a^)x- (4) 

Let us seek multi-phase solutions in the form a = a{R^ , R^), b = b{R^ , R"'), p = 
p{R^ , R"") where the 'phases' R^{x,y,t) satisfy the equations (2). The substitution of 
this ansatz into (4) implies the relations 

dib = iJdia, dip = X^dia, A* = a + (iJf. 

Calculating the compatibility conditions didjb = djdib, didjp = djdip, and substituting 
= a + {iJi^Y into the commutativity conditions (3), one obtains the Gibbons-Tsarev 
system for a{R) and fJ.\R), 

r, i dja ^ ^ ^ diadna 



i\2 ' 
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i ^ ii which was first derived in [26] in the theory of hydrodynamic reductions of Benney's 
moment equations. It is remarkable that the Gibbons-Tsarev system is in involution, and 
its general solution depends, modulo reparametrizations i?* on n arbitrary 

functions of one variable. Thus, the dKP equation possesses infinitely many n-component 
reductions parametrized by n arbitrary functions of one variable. We point out that the 
compatibility conditions dkdjH^ = djdkH^ and didjdka = didkdja involve triples of indices 
^ 7^ J 7^ ^ only. Thus, for = 2 the Gibbons-Tsarev system is automatically consistent, 
while its consistency for n = 3 implies the consistency for arbitrary n. Based on this 
example, we give the following 

Definition An equation of the form (1) is said to he integrable if, for any n, it pos- 
sesses infinitely many n-component hydrodynamic reductions parametrized by n arbitrary 
functions of one variable. 

We have verified that all examples presented above are indeed integrable in this sense, 
with the exception of the equations Hess u = 1 and Hess u = Am, which do not pass 
the test (see Sect. 3.5). In Sect. 2 we derive the integrability conditions (as a system of 
third order differential relations for the function F in (1)), and prove our first main result 
(Theorem 1 of Sect. 2): 

• The moduli space of integrable equations of the dispersionless Hirota 
type is 21-dimensional. 

The class of equations (1) is form-invariant under the action of the contact group 5*^(6) 
generated by linear symplectic transformations of the variables x,y,t,Ux,Uy,Ut. These 
transformations map integrable equations to integrable. Taking into account that dim 
Sp{G) = 21, our second main result (Theorem 5 of Sect. 6) reads as follows: 

• The action of the equivalence group Sp{Q) on the moduli space of inte- 
grable equations has an open orbit. 

This result is, in a sense, surprising: it estabhshes the existence of a unique 'master- 
equation' which generates all other integrable examples via various (singular) limits. The 
existence of an open orbit does not mean, of course, that all integrable equations of the 
form (1) are S'p(6)-equivalent. The structure of the orbit space of this action is quite 
complicated, for instance, all examples listed above belong to various singular orbits of 
lower dimension, and are not S'p(6)-equivalent. 

From differential-geometric point of view, the equation (1) can be viewed as defining 
a 5-dimensional hypersurface in the Lagrangian Grassmannian A which can (locally) 
be identified with the space of 3 x 3 symmetric matrices Uij. In Sect. 4 we demonstrate 
that each hypersurface C A inherits a special structure, namely, a field of rational 
normal curves 7 of degree four specified in the tangent bundle TM^. Bisecant surfaces 

C are defined as two-dimensional submanifolds whose projectivised tangent spaces 
are bisecant lines of 7. Similarly, trisecant 3-folds are three-dimensional submanifolds 

C whose projectivised tangent spaces are trisecant planes of 7. Our third main 
result is the following geometric characterization of integrable equations (Theorem 3 of 
Sect. 4): 
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• Bisecant surfaces and trisecant 3-folds are geometric images of the two- 
and three-component hydrodynamic reductions, respectively. In partic- 
ular, the equation is integrable if and only if the corresponding hyper- 
surface possesses infinitely many trisecant 3-folds. 

The paper is organized as follows. 

In Sect. 2 we derive the integrability conditions by transforming (1) to a quasilin- 
ear form, and calculating the compatibility conditions of the corresponding generalized 
Gibbons-Tsarev system which governs n-component hydrodynamic reductions. The re- 
sult of this calculation is a complicated set of the integrability conditions, which are in 
involution. This establishes the 21-dimensionality of the moduli space. 

Further examples of integrable equations and partial classification results are provided 
in Sect 3, based on the integrability conditions obtained in Sect. 2. In particular, we 
analyze the integrability of the symplectic Monge- Ampere equations. 

Differential-geometric aspects of integrable equations of the form (1) are discussed in 
Sect. 4. We begin by introducing a flat generalized conformal structure on the Lagrangian 
Grassmannian A, whose conformal automorphism group is isomorphic to 5*^(6). Given a 
hypersurface C A, we describe geometric structures which are naturally induced on 
M^, and provide a geometric characterization of the integrability conditions. 

In Sect. 5 we briefly discuss multi-dimensional equations of the form (1). 

The action of the equivalence group 5*^(6) on the moduli space of integrable equations 
is discussed in Sect. 6. 

2 Derivation of the integrability conditions: proof of 
Theorem 1 

In this section we derive the integrability conditions which follow from the requirement 
of the existence of n-component reductions, and prove the following 

Theorem 1 The moduli space of integrable equations of the dispersionless Hirota type is 
21- dimensional. 

Proof: 

To apply the method of hydrodynamic reductions we flrst rewrite the equation (1) in 
the evolutionary form, 

"^tt — fi^xxy Uxy, Uyy, Uxt, Uyt) , (S) 

and introduce the notation 

Uxx = a, Uxy = b, Uyy = c, Uxt = p, Uyt = q, Uu = f{a, b, c,p, q). 
This provides an equivalent quasilinear representation of (5), 

ay = bx, at = Px, by = Cx, h = Py = qx, Ct = qy, 

(6) 

Py = Qx, Pt = fia,b,c,p,q)x, qt = f{a,b,c,p,q)y. 



5 



Looking for solutions in the form a = a{B}, -R"), h = b{R^, i?"), c = c{R^, K^), p = 
p{R^, K^), q = q{R^ , R"'), where the Riemann invariants R^ satisfy the equations 
(2), and substituting this ansatz into (6), we obtain 

dib = f/diU, dip = X^dia, dig = n^^dia, diC = {i/^dia, (7) 

di = d/djit, along with the dispersion relation 

D(A„ /iO = fa + ftfi' + Ufi^ + fpX' + /,A>^ - (A^^ ^ 0; (8) 

in what follows we assume that the dispersion relation (8) defines an irreducible conic in 
the (A, /i)-plane. This is equivalent to the requirement that the expression 

^ fb ~^ fbfpfg fafq fcfp ^fafc 

does not vanish. Calculating the consistency conditions for (7) we obtain 



iOjtt \j _\i OiQj + -j-j—j-j-OjCl, 



djXd-ia + diX^dja = 0. 



(9) 



Applying the operator dj to the dispersion relation (8) and using (7) and (3), we obtain 
(9jA* and djfi^ in the form 

djX = (A* - X^)Bijdja, djfi' = (/i* - fi^)Bijdja (10) 

where Bij are rational expressions in A*, A-', /i', fi^ whose coefficients depend on the second 
order partial derivatives of the function f{a,b,c,p,q). Explicitly, one has 

where 

+ fbcfi'fi'ifi' + /i^) + fbpiX'fi^ + AV') + /6,/i>^(A' + A^) 

(11) 

+ + fcp{y{fl^f + X^ifi')') + /c,yU>^(A>^- + AVO 

+/ppA^A^' + /p,A^A^(/i^ + fi^) + /,,AUV>^' 

and 

Aj = -2yy + 2/„ + + + 2/,f,>J + f,(y + \') + /,(a>' + a>') 

^,d(^,^±^, (12) 



the last equality holding modulo the dispersion relation (8). Substituting (10) into (9)i 
one obtains 

didjtt = —2Bijdiadja; (13) 

the second condition (9)2 holds identically by virtue of the symmetry Bij = Bji. The 
compatibility conditions dkdjX^ = djO^X^, dkdjfi^ = djdkfi^ and dkdjdia = djdkdia are 
equivalent to the relations 

dkBij = {BijBkj + BijBik — BkjBik)dka, (14) 

which must be satisfied identically by virtue of (7), (8) and (10). 

Remark. Notice that each compatibility condition involves three distinct indices only. 
This observation immediately implies that 

(i) any equation of the form (1) possesses infinitely many two-component reductions 
parametrized by two arbitrary functions of one variable. Indeed, these reductions are 
governed by the equations (10) and (13) where i,j = 1,2 and A*,/i* satisfy the disper- 
sion relation (8). These equations are automatically consistent, and the general solution 
depends, modulo reparametrizations i?* — ^ ip'^{R^), on two arbitrary functions of one vari- 
able. Therefore, the existence of two-component reductions is a common phenomenon 
which is not related to the integrability. 

(ii) the existence of three-component reductions implies the existence of n-component 
reductions for arbitrary n. Thus, one can define the integrability as the existence of 
infinitely many three-component reductions parametrized by three arbitrary functions of 
one variable. This property is reminiscent of the well-known three-soliton condition in 
the theory of integrable systems. 

In order to simplify the derivation of the integrability conditions we first rewrite (14) 

as 

dkNij = Nij (^-^dkDij + Bkjdka + Bikdkoj - DijBkjBiAa. (15) 

Here, the third order derivatives of /(a, 6, c, p, q) are present only in the l.h.s. term dkNij 
(see the form of Nij and Dij above). Further reduction of the complexity of the expression 
in the r.h.s. is achieved by representing 1/Dij in the form 

= Uij = [2fa + (/il + H2)fb + 2/il/i2/c - (Al + A2)/p - (Ai/i2 + A2/il)/g 
+ fp + (/^l + f^2)fpfg + /il/i2/,' + 2AiA2]/((/il - fi2?A) 

(which holds identically modulo the dispersion relation (8)), and the subsequent substi- 
tution Bst = Nst/Dst = NstUst- The denominators of the r.h.s. terms in (15) cancel out, 
producing a polynomial in \\ , fi\ fi^ , fi^ with coefficients depending on the deriva- 
tives of the density f{a,b,c,p,q). This was the crucial simplification of the calculation: 
the starting expression for the r.h.s. of (15) has more than 450.000 terms with different 
denominators; after properly organized cancellations we get a polynomial expression with 
less than 6.000 terms. Using the dispersion relation (8), we simplify this polynomial by 
substituting the powers of (A*)^, (A-')*, (A'^)'*, s > 2, arriving at a polynomial of degree 
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one in each of A*, A-', A'^, and degree two in /i's. Equating similar coefficients in both sides 
of (15), we arrive at a set of 35 equations for the derivatives of the function /(a, b, c,p, q), 
which are hnear in the third order derivatives. Solving this linear system we obtain the 
closed form expressions for all third order derivatives of /(a, b, c,p, q) in terms of its ffist 
and second order derivatives, which we represent symbolically in the form 

d'f = Ridf,d'fy, (16) 

here R denotes complicated rational expressions in the ffist and second order derivatives 
of / (it is clear from the above that the function / itself will not enter these expressions 
explicitly). A straigtforward computer calculation shows that the overdetermined system 
(16) is in involution. Thus, the moduli space of integrable systems of the type (1) is 
21-dimensional: one can arbitrarily specify the values of /, df and d'^f at any fixed point. 
This amounts to 1 + 5 + 15 = 21 arbitrary constants. This finishes the proof of Theorem 
1. 

The right hand sides of (16) are not presented here because of their complexity (see 
Sect. 3 for a discussion of particular cases where these formulae become less cumbersome, 
and are presented explicitly). The integrability conditions (16) provide a straightforward 
computer test of integrability for any equation from the class under consideration, and 
allow one to obtain classification results. These conditions also allow computer checks 
of various properties of differential-geometric objects naturally associated with equations 
from the class (1). 

Remark. In two dimensions, any second order PDE of the form F{uxx, u^y, Uyy) = 
is automatically integrable. Indeed, introducing the parametrization 
6, Uyy = f{a,b), so that F{a,b, f{a,b)) = 0, one obtains a two-component quasilinear 
system 

ay = h, by = f{a,b)x; 

any system of this type linearises under a hodograph transformation which interchanges 
dependent and independent variables. This simple trick, however, does not work in more 
than two dimensions. 



3 Examples and classification results 

In this section we discuss examples (both known and new) and partial classification results 
of integrable equations of the form (1) which are obtained using the integrability conditions 
(16) derived in Sect. 2. 

3.1 Integrable equations of the form Uu = f{uxx,Uyy) 
In this case the integrability conditions (16) simplify to 

r f 1 •^'^'^ I I / f I ^'^'^ _L ^'"''^ 

Jaaa Jaa I j, ~r ~^ I , Jaac Jaa I j, ~r 



c Ja J \ Jc Ja 

r r I fee fac \ n f ( ^'^'^ i '^'^'^ | / / ( f ^"^ 

face fee \ ~7 I 7~ I ; feee fee I ~7 I 7~ J ; faafce {fc 

fc fa / \ fe fc 



ae J 1 
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recall that, according to the notation introduced in Sect. 2, This 
system is in involution and its general solution depends on 5 integration constants. To 
solve it explicitly one notices that the first two equations imply = const. Similarly, 

Ja Jc 

fee 
fafc 

following general solution. 



the next two equations imply = const. Further elementary integration leads to the 



/(a, c) = s ln(me'"' + ne"") + const. 
This corresponds to equations of the form 

^gM«.. + ne'^'^yy + fee''"" = 0, 
as well as degenerations thereof, 

where the coefficients are arbitrary constants. In fact, all these coefficients can be elimi- 
nated by appropriate complex rescalings, leading to the two essentially different examples, 

g«.. _^ ^Uyy ^ ^Utt g^j^^ _^ ^ g«tt_ 

The first equation is apparently new, while the second is the Boyer-Finley equation. 

3.2 Integrable equations of the form Uu = f{uxx,Uxy,Uyy) 

This is a generalization of the previous case. Thus, we assume fb 7^ 0. The integrability 
conditions (16) take the form 

faafbb fab i faafcc fac 5 fbbfcc fbc 1 

faafbc fab fact fab fee faefbei fab fbc facfbbj 

/— J au, yjujuu —jcjaa —jajae; p ^fab {fbfab ^fcfaa ^fafac) 

aaa 7"2 ~r ~r ' Jaab „ 2 



fb ^fafe 

'2 fab {fbfbb — '^fefab — '^fafbc) 
Jaac r ' J abb „ 2 



'^faa {fbfab 


"^fcfaa 


~ '^fafacj 




f? 


^fafc 




Vac 


■ {fbfab ~ 


'^fcfaa ~ 


'^fafacj 




fb' 


"^fafc 




'2 fac 


{fbfbc — 


'^fcfac 


'^fafccj 




fb' 


'^fafc 




'^fbb 


{fbfbb — 


'^fcfab 


'^fafbc) 




fb' 


^fafc 




'2 fbc 


{fbfbc — 


'^fcfac 


'^fafcc) 



/~ J ac \j uj uc —jcjac —jajcc/ p 
acc r"2 r ' J' 



I abc 



fb ^fafc 

'2 fab {fbfbc '^fcfac '^fafa 
2 



fb '^fafc 



n ^ J UU \j UJ UU —jcjau —jajuc; j. "2 fbc {fbfbb — '^fcfab — '^fafbc) 

Jbbb c 2 J c c ' Jbbc „ 2 



fb ^fafc 
"2 fee {fbfbc '^fcfac '^fafa 



r ^ J uc \ j UJ uc ^ J cj ac ^ j aj cc j n 

^ fb-^fafc ' ~ fb' - ^fafc 
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Notice that the first six equations imply that the Hessian matrix of the function /(a, b, c) 
has rank one, so that one can set fa = m{fi)), fc = n{f}y). The substitution into the 
remaining equations imphes either m" = n" = 0, or fbb = 0. Further elementary analysis 
leads, up to linear transformations of x and y, to the following canonical forms (here we 
only list those representatives which contain a nontrivial dependence on Uxy)'- 

(a) Utt = CiUxx + P'Uyy + fiur^y), here ip satisfies a third order ODE ip"'{{ip'Y — iaP) = 
2ip'[ip"y. The integration leads to the three canonical forms, 

Utt = au^x + /3uyy + In cosh -fu^y, 

Utt = auxx + l^^Uxy, or utt = \nuxy, 

the last two cases corresponding to /5 = and a = P = 0, respectively. The first two 
examples can be viewed as generalizations of the Boyer-Finley equation. 

(b) Utt = u^y + f3{u^x - Uyy) + if{u^^ + Uyy), lieie (/?'"(4(y?')^ - 1 - 4/?2) = 8if'{if")'^. This 
leads to the equation 



Utt = Uxy + f3{UxX - Uyy) H hi cosh. '-i{UxX + Uyy). 



1 + 4/32 
27 



(c) Utt = Uxy + Puxx + f{uyy), here ip satisfies a third order ODE ip"'{A[3ip' — 1) = 4,[3{ip")'^ . 
This results in the equation 

Utt = Uxy + (3uxx + ^^ra + «e^"''^ 
whose degeneration, corresponding to /3 = 0, is the dKP equation Utt = u^y + u^y. 

3.3 Integrable equations of the form Uxy = f{uxt, Uyt) 

Formally, equations from this class are not of the form (5), however, they can readily be 
made 'evolutionary' by an appropriate linear change of the independent variables. The 
resulting set of integrability conditions looks as follows: 

f — f ( hn I Ipp_\ f _ f I Isg I fpg 
Jppp — Jpp \^ J. ~^ j J ^ JPPi ~ JPP y f f 

f — f ( im -L Iee.] f — f ( isi _i_ Im 

Jpqq Jqq \^ j. j J ^ Jill J 11 \^ j f 

This system is in involution, and its general solution depends on 6 integration constants. 
To solve it explicitly one notices that the first two equations imply = const. Similarly, 

p q 

the last two equations imply -j^ = const. Further elementary integration gives, under 
the assumption that both fpp and fgq are nonzero, the general solution 

f{p,(l) = - In 



K \ re'^p+'^i 
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which leads to the dipersionless Hirota-type equation for the BKP hierarchy [7], 
Up to complex rescalings, it can be transformed to 

Degenerations of (17), corresponding to fpp = or fpp = fgg = 0, result in 

Uxy = u^t + e"^*, = u-j^t tan(nj^t) 

and 

respectively. 

3.4 Integrable equations of the form Uu = f{uxx,Uxt,Uxy) 

Equations of this type arise in the theory of integrable hydrodynamic chains satisfying 
the additional' Egorov' property. Here we reproduce the classification result from [41] (see 
also [11], [25]). The integrability conditions (16) take the form 



fbbb — 


2Ai 

fb ' 


fabb 


'^fabfbb 
fb ' 


fpbb — 


'^fpbfbb 
fb 


faab 


V!b 
fb ' 


fapb 


'^fabfpb 
fb ' 


fppb ~ 


Vlb 

fb ' 



2 

fppp — "72 {fpfpb ~^ fpbifbfpp + 2/at,) fbb{fpfpp + '^fap)^ ; 
Jb 

2 2 

fapp '72 i^fafpb ~^ fabifbfpp ~l" fab) fbbifafpp ~l" faaS) ) (-'-^) 

h 

2 

faap 'T2 ifbbifpfaa '^fafap) fabifpfab '^fbfap) fpbifbfaa fa fab)) y 

h 

faaa = ~P2 {{fa + fp) fab ~^ fafpb ~^ fb ifap ~ faafpp) ~ fppfbbfa 
Jb 

~\~fabfb{faa ~\~ '^{fafpp fpfap)) + '^fpb{fp{fbfaa fafab) fafbfap) 
fbbiifa ~l" fp)faa "^fafpfap)^ ; 

this system is in involution and its general solution depends on 10 arbitrary constants. 
The integration of these equations leads to the four essentially different canonical 
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forms, 



Ucr.cr. V Uxx ^^xx J ^xx '^^xx 



Utt = — + ■^V{Uxx)ult, 

Uxt 6 



Utt 



Inu^y - \n9i {u^t, u^x) ~\j v{'r)dT, 

see [41]. Here A, B, C are arbitrary constants, 77 is a solution to the Chazy equation [16], 

7]"' + 2r]r]" = 3(77')^, 

and 61 is a theta-function in the variable Uxt whose dependence on Uxx is governed by the 
Chazy equation. 

3.5 Symplectic Monge- Ampere equations 

Let us consider a function u{x^, ...x^) of k independent variables, and introduce the x 
Hessian matrix U = [uij] of its second order partial derivatives. The symplectic Monge- 
Ampere equation is a PDE of the form 

Mfc + Mk-i + ... + Ml + Mo = 

where Mi is a constant-coefficient linear combination of all / x / minors of the matrix 
U, < I < k. Thus, Mfc = det U = Hess u, Mq is a constant, etc. Equivalently, these 
PDEs can be obtained by equating to zero a constant-coefficient /c-form in the 2k variables 
x\ Ui. This class of equations is invariant under the natural action of the symplectic group 
Sp(2k). In the case k = 2 one obtains a standard Monge- Ampere equation, 

U11U22 - M12 + ctuu + l3ui2 + 7M22 + 5 = 0, (19) 

which can be interpreted as the equation of a 'sphere' corresponding to the pseudo- 
Euclidean metric duudu22 — du\2- Monge-Ampere equations (19) can be characterized as 
the only equations of the form -F(tiii, iii2, ^22) = which are linearizable by a transfor- 
mation from the equivalence group 5*^(4). 

The case /c = 3 is also understood completely: one can show that, for /c = 3, any 
symplectic Monge-Ampere equation is either linearizable (in this case it is automatically 
integrable), or S'p(6)-equivalent to either of the two essentially different canonical forms 
[37, 4], 

Hess M = 1, Hess u = uu + U22 + Mss- (20) 

Based on the integrability conditions (16), we have verified directly that both PDEs (20) 
are not integrable by the method of hydrodynamic reductions. Thus, a 3-dimensional 
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symplectic Monge-Ampere equation is integrable if and only if it is linearizable (this is 
no longer true in more than three dimensions - see Sect. 5). The linearizability condition 
constitutes a single relation among the coefficients of the equation: for a Monge-Ampere 
equation of the form 



e det 



Mil Ui2 Ui3 
Ui2 U22 U23 
Ul3 U23 M33 



+ hi{u22U33 - wis) + h2{UuU33 



U 



13, 



+ h3{uiiU22 - U 



2 1 

12) 



+9l{unU23 - U12U13) + g2{u22Ul3 - U12U23) + g2,{u33Ui2 - U13U23) 



(21) 



+ S2U22 + S3U33 + T1U23 + T2U13 + T3U12 + u = 0, 
the linearizability condition defines a quartic hypersurface in the space of coefficients, 

hfsl + hlsl + hlsl + gls2S3 + glsiS3 + glsiS2 

-2(/ii/l2SiS2 + hih3SiS3 + h2h3S2S3) + AeSiS2S3 + Al/hih2h3 

+enT2T3 - vgig2g3 - e^'^^ - i^(5'?^i + 5-2^2 + 5'|^3) 
-(fi'iTi + ^2^-2 + g'iT3 + 2eij){hisi + h2S2 + h3S3 - ev) (22) 

+2{gihiSiTi + 5-2/^252^2 + g3h3S3T3) + /l2/i3 + r| /il/i3 + r| /il/i2 
-e(rf Si + r|s2 + ^-fsg) + SiTig2g3 + S2T2gig3 + SsTg^i^fa 
-{gihlT2T3 + g2h2TiT3 + g3h3TiT2) = 0. 

Notice that, if e 7^ 0, one can always eliminate second order minors in (21) by adding to 
u an appropriate quadratic form. In this case the equation takes the form 



det 



un U12 Ui3 
U21 U22 U23 
U31 U32 U33 



+ SiMii + S2M22 + S3U33 + T1U23 + r2Ml3 + r3Mi2 + Z/ = 0, 



while the linearizability condition simplifies to 

4S1S2S3 + i/^ + T1T2T3 - SitI - S2T2 - 53X3 = 0. 

Remark 1. The condition (22) can be given an invariant formulation [37, 4]. Let Q = 
dui A dx^ be the standard symplectic form, and v = vq the associated canonical bivector. 
Let u! be the effective 3-form corresponding to the equation (21): recall that u is effective 
if a; A r2 = 0. Let us define the quadratic form 



(23) 
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where ix^ is the inner product of uj with a vector X, and _L is the operator of inner 
multiphcation by the canonical bivector. The hnearizabihty condition (22) is equivalent 
to the requirement that the quadratic form (23) is degenerate. 

Remark 2. The equation (20) i arises in the theory of improper affine spheres, while the 
equation (20)2 describes special Lagrangian 3-folds in [13, 32]. Their contact-equivalent 
forms are 

2 

and 

2 I 2 I 2 1 

UxxUyy '^xy ' '^XX^tt xt ' '^ttUyy Itjy — i, 

respectively (both non-integrable). 

Remark 3. Equations of the form (21) and the hnearizabihty condition (22) have a 
clear projective- geometric interpretation. Let us consider the Lagrangian Grassmannian 
A, which can be (locally) identified with the space of 3 x 3 symmetric matrices U (see 
Sect. 4 for more details). The minors of U define the Veronese embedding of A into 
the projective space P^^ (we will identify A with its projective embedding). Thus, sym- 
plectic Monge-Ampere equations correspond to hyperplane sections of A. Linearizable 
equations correspond to hyperplanes which are tangential to A. Thus, the hnearizabihty 
condition (22) coincides with the equation of the dual variety to A, which is known to be 
a hypersurface of degree four in (p^^-j* (^^^^ j-gg^ gg^ ^gj ^ general theory behind this 
example) . 

The results of this section answer in the negative the question formulated by Joyce 
in [32]: do special Lagrangian m-folds in C™ for m > 3 constitute some kind of higher- 
dimensional integrable systeml 

4 Differential-geometric aspects of integrable equa- 
tions of the dispersionless Hirota type 

In this section we adopt a geometric point of view, and consider (1) as the equation 
difining a 5-dimensional hypersurface in the 6-dimensional space with coordinates Uij 
(which, as we will see shortly, is naturally identified with the Lagrangian Grassmannian 
A). Our aim is to reformulate the integrability conditions in the intrinsic geometric terms. 

In Sect. 4.1 we discuss differential geometry of the Lagrangian Grassmannian A. We in- 
troduce the main object needed to develop differential geometry of hypersurfaces, namely, 
a canonically defined symmetric cubic form which is invariant under the natural action of 
the symplectic group Sp{Q). This form (to be precise, its conformal class), plays the role 
of an ambient flat conformal structure. 

In Sect. 4.2 we develop differential geometry of hypersurfaces C A. The induced 
cubic form gives rise to a field of rational normal curves of degree four defined in the 
projectivization of the tangent bundle TM^. Thus, a hypersurface of the Lagragian 
Grassmannian carries an intrinsic irreducible GL(2)-structure. 



14 



In Sect. 4.3 we show that two-component and three-component hydrodynamic reductions 
of the equation (1) correspond to bisecant and trisecant submanifolds of M^, respectively. 
This allows one to reformulate the integrability as the existence of sufficiently many 
trisecant submanifolds parametrized by three arbitrary functions of one variable. 

Some useful computational formulas are provided in Sect. 4.4. 

General aspects of irreducible GL(2)-structures in dimension 5 are discussed in Sect. 4.5. 
We demonstrate that the requirement of the existence of bisecant surfaces imposes strong 
constraints on the torsion of the corresponding GL(2)-structure. 

4.1 Sp{6) as a symmetry group of the problem: differential ge- 
ometry of the Lagrangian Grassmannian 

Let us consider a 6-dimensional symplectic space with canonical coordinates x = (x^, x^, x^Y 
and p = {pi,P2,P3y, viewed as 3-component column vectors. Lagrangian submani- 
folds can be parametrized in terms of a generating function u{x,y,t): x = {x,y,ty and 
p = [ux, Uy, utY. Lagrangian planes are defined by the equation cip = f/rfx where [/ is a 
3x3 symmetric matrix (the Hessian matrix of u). Thus, the Lagrangian Grassmannian 
A is 6-dimensional, and can (locally) be identified with the space of 3 x 3 symmetric ma- 
trices. The equation (1) defines a 5-dimensional hypersurface C A; the corresponding 
solutions u{x, y, t) can be interpreted as Lagrangian submanifolds whose Gaussian images 
belong to the hypersurface M^. 

The action of the linear symplectic group Sp{6), 

f dp\^f A B\f dp\ 
\dSt J \C D J \ J' 

naturally extends to A: 

U={AU + B) {CU + D)-\ (24) 

Here A, B, C, D are 3 x 3 matrices such that A^C = C^A, B^D = D^B, A^D - C^B = 
id] notice that the extended action is no longer linear. The transformation law (24) 
suggests that the action of 5*^(6) preserves the class of equations (1), indeed, second 
order derivatives transform through second order derivatives only. Moreover, since any 
changes of variables obviously preserve the integrability, the group Sp{Q) can be viewed 
as a natural equivalence group of the problem: it maps integrable equations to integrable. 
Thus, 5*^(6) is a point symmetry group of the integrability conditions (16) derived in Sect. 
2. The classification of integrable equations of the form (1) has to be performed modulo 
this equivalence: two S'p(6)-related equations should be regarded as 'the same'. 

Geometrically, our problem is reduced to the classification of hypersurfaces of the 
Lagrangian Grassmannian A (satisfying certain 'integrability' conditions to be specified 
later, see Sect. 4.3), up to the action of 5*^(6). Our next goal is to clarify which differential- 
geometric structures are induced on hypersurfaces of A. To do so one first needs to 
introduce differential-geometric objects which are naturally defined on the Lagrangian 
Grassmannian, and are invariant under the action of 5*^(6). 
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The condition of nontrivial intersection of two infinitesimally close Lagrangian planes 
defined by symmetric matrices U and U + dU, is det dU = 0. This condition is manifestly 
invariant under the action of Sp{6) as specified by (24). Thus, each tangent space to the 
Lagrangian Grassmannian A is equipped with a cubic cone C defined by the equation 
det dU = 0. The projectivisation of this cone is known as the 'cubic symmetroid'. The 
singular locus V C C of the cubic cone is specified by the condition rank dU =1. Its 
projectivization is known as the Veronese variety, which is a non-singular algebraic surface 
of degree four. 

The identity 

det{CU + D) 

which readily follows from (24), implies that the conformal class of the cubic form det dU 
is invariant under the action of 5*^(6). The converse is also true: 

Theorem 2 The group of conformal automorphisms of the symmetric cubic form det dU 
is isomorphic to Sp{Q). 

The proof consists of a direct calculation of conformal automorphisms of the cubic form 
det dU. The corresponding infinitesimal generators are presented in Sect. 6. In fact. 
Theorem 2 can be traced back to the old paper by Cartan [15], Theorem XX. We refer to 
[5, 27] for further generalizations of this Liouville-type result. 

Remark. The theorem is true in any dimension. Thus, for the Grassmannian of 
Lagrangian 2-planes in a 4-dimensional symplectic space (naturally identified with 2x2 
symmetric matrices U), the expression det (it/ defines the Lorentzian metric duiidu22 — 
du{2 whose conformal automorphisms generate the group S'0(3,2). This establishes a 
well-known isomorphism between Spi^A) and 5*0(3,2). 



4.2 Geometric structures on a hypersurface of the Lagrangian 
Grassmannian 

Let be a hypersurface of the Lagrangian Grassmannian A. Taking a point s G 
and intersecting the tangent space TgM^ with the singular locus V of the cubic cone C 
in T5A, one obtains, after a projectivisation, a rational normal curve of degree four (it 
is well-known that a generic hyperplane section of the Veronese variety \^ is a rational 
normal curve). This curve can also be interpreted as a set of matrices of rank one in 
the tangent space TgM^ (recall that T^A is identified with a space of 3 x 3 symmetric 
matrices). Thus, the projectivised tangent bundle of the hypersurface is equipped 
with a field of rational normal curves. Since the group of conformal automorphisms of 
the rational normal curve is isomorphic to GL{2), this specifies on an irreducible 
GL(2)-structure. Alternatively, one can say that each tangent space to is canonically 
identified with a five- dimensional space of binary quartics; in this picture the rational 
normal curve corresponds to quartics with a quadruple root. Similar structures were 
discussed in [19, 3] under the name 'paraconformal'. 
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Remark. Another natural source of irreducible G'L(2)-structures in dimension five 
is provided by 5-th order ordinary differential equations with the vanishing Doubrov 
invariants [17, 19]. Here GL(2)-structures arise on the moduli spaces of solutions to the 
corresponding ODEs. The relation of these structures to equations of the dispersionless 
Hirota type is yet unclear. 

Irreducible GL(2)-structures in dimension five possess a number of specific properties, 
allowing one to equip with more 'familiar' objects such as the conformal cubic form 
and the conformal metric. These are defined as follows: 

The conformal cubic form Ca/s is a restriction to of the cubic form C defined 
on A. The knowledge of this cubic form allows one to reconstruct the GL(2)-structure, 
indeed, the singular locus of the cubic C^s is the rational normal curve introduced above. 
Conversely, given a rational normal curve, its bisecant variety is the cubic Cjvfs = [29], 
p. 119. This specifies Cm^ uniquely up to a conformal factor. 

The conformal metric Qm^ is defined by the unique quadratic cone in TgM^ containing 
the tangent variety of the rational normal curve, see [29], p. 119, for an alternative defini- 
tion. The tangent variety of the rational normal curve is a complete intersection of Qm^ 
and C^s. 

One can show that the conformal metric Qm^ and the cubic form C^/s satisfy a system 
of remarkable relations which have appeared recently in [6] (notice that the signature of 
the quadratic form in [6] is Euclidean, while in our case the metric is pseudo-Euclidean 
of the signature (3,2)). Let Qij and Cijk be the components of Qms and Cm^ in any 
coordinate system on (we point out that these objects are defined up to conformal 
factors). One can show that the conformal factors can be specified in such a way that the 
following relations are satisfied: 

CijkQ^-^ = (apolarity condition), 

(25) 

CjkrQ Clns ~\~ CljrQ C^ns ~l~ C^lrQ Cjns QjkQln ~l~ QljQkn ~l~ QklQjm 

these relations still allow simultaneous conformal reparametrizations of the form C 
(f^C, Q — s> ip'^Q. The explicit coordinate formulas for Qm^ and Cms are provided in Sect. 
4.4. 

The normalised conformal metric Qm^ said the cubic form C^s satisfying the equations 
(25) carry important geometric information about the hypersurface M^. Their transfor- 
mation laws suggest that one should consider the ratio a = {Cm^YHQm^Y which can 
be viewed as an analog of the Fubini projective element. By construction, this object is 
manifestly S'p(6)-invariant. We can formulate the following 

Conjecture. A generic five- dimensional hypersurface of the Lagrangian Grassmanian 
A is uniquely defined, modulo Sp{Q)- equivalence, by its 'symplectic element' a . 
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4.3 Geometric interpretation of the integrability conditions: bise- 
cant surfaces and trisecant 3-folds 

As explained in Sect. 4.2, each projectivised tangent space TgM^ carries a rational normal 
curve 7 of degree four. Among the most natural geometric objects associated with a 
rational normal curve are its bisecant lines (two-parameter family) and trisecant planes 
(three-parameter family) . 

Definition. A bisecant surface is a two-dimensional submanifold C whose projec- 
tivised tangent planes are bisecant lines. Similarly, a trisecant 3-fold is a three-dimensional 
submanifold C whose projectivised tangent spaces are trisecant planes. To be more 
precise, we will consider holonomic trisecant 3-folds which can be defined as follows. No- 
tice first that each tangent space T^S^ carries three distinguished directions, namely, those 
corresponding to the three points of intersection of PTsT? with 7. These directions define 
a net on T?, which we will call the characteristic net. We require the characteristic net 
to be holonomic (that is, a coordinate net). A similar net exists on any bisecant surface, 
however, the requirement of holonomicity is superfluous: in two dimensions, any net is 
automatically holonomic. 

Theorem 3 Bisecant surfaces and holonomic trisecant 3-folds of a hypersurface cor- 
respond to two- and three- component hydrodynamic reductions of the associated equation 
of the dispersionless Hirota type. Moreover, 

(i) Each hypersurface possesses infinitely many bisecant surfaces parametrized by two 
arbitrary functions of one variable; 

(ii) A hypersuface corresponds to an integrable equation if and only if it possesses 
infinitely many holonomic trisecant 3-folds parametrized by three arbitrary functions of one 
variable. Thus, the existence of holonomic trisecant 3-folds is a geometric interpretation 
of the integrability property. 



We follow the notation of Sect. 2. Consider an equation of the dispersionless Hi- 
rota type represented in the form (5), and take its n-component reduction specified 



by a = a{R\...,R'), b = b{R\ R""), c = c{R\ R""), p = p(i?\ i?"), q = 



q{R^ , R^), f = f{R^,...R^). The image of this solution in the Lagrangian Grass- 
mannian is a submanifold E"' C (it would be sufficient for our purposes to restrict to 
n = 2, 3), represented by a symmetric matrix 



parametrized by Using (7) and (8), one obtains the following expression for 

the derivative of U with respect to R\ 



Proof: 




(26) 




(27) 
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thus, Ui is a matrix of rank one, so that the projectivization of diU belongs to 7, and the 
coordinates i?* provide a characteristic net on S". For n = 2 we have a two-dimensional 
surface parametrized by R^, B?. Since both diU and have rank one, the surface T? 
is bisecant. As follows from Sect. 2, any equation of the dispersionless Hirota type (not 
necessarily integrable) possesses infinitely many two-component reductions parametrized 
by two arbitrary functions of one variable. Thus, any hypersurface of the Lagrangian 
Grassmannian possesses an infinity of bisecant surfaces. This establishes the first part (i) 
of the theorem. 

Similarly, any three-component reduction corresponds to a holonomic trisecant 3-fold 
C parametrized by R^,R^,R'^. Since an integrable equation possesses (by defini- 
tion) infinitely many three-component reductions parametrized by three arbitrary func- 
tions of one variable, the corresponding hypersurface possesses an infinity of holonomic 
trisecant 3- folds. This establishes the first part of (ii). To finish the proof one needs to 
show that, conversely, bisecant surfaces (holonomic trisecant 3-folds) of corresponds 
to two-component (three-component) reductions of the associated equation. This can 
be demonstrated as follows. Let be a bisecant surface represented in the form (26), 
referred to its characteristic net R^,R^. Thus, a,b,c,p,q, f are functions of R^,R^ such 
that the rank of diU equals one, so that one can introduce the parametrization (27). The 
compatibility conditions for the equations dib = fi'^dia and dip = X^dia imply 

did^a = rditt H — rd^a and &9,a = — r-fta + — r-dia, (28) 

respectively. Using (28), the consistency conditions for diC = {jj^Ydia and dif = {X^ydia 
simplify to djfi^dia+difi^dja = and djX^dia+diX^dja = 0, compare with (9). Substituting 
these relations into (28) one obtains 



didjtt = 2 — — rdia and didja = 2—^ — r-dia, 

1^^ - jJ'' AJ - A* 



which implies the commutativity condition (3). Finally, the last compatibility condition, 
diQ = fi^X^diQ, holds identically. The dispersion relation (8), which connects /i* and 
A*, follows from the relation d^f = {X^Ydia after the substitution / = f{a,b,c,p,q). 
Ultimately, we have recovered all equations defining hydrodynamic reductions, see Sect. 
2. Thus, bisecant surfaces correspond to two-component reductions (for trisecant 3-folds, 
considerations are exactly the same). 

We have characterized the integrability as the existence of 'sufficiently many' special 
submanifolds. This is very reminiscent of the characterization of self-dual 4-manifolds in 
terms of the existence of a-surfaces [1], or a similar characterization of semi- integrable 
almost Grassmann structures [2]. 

4.4 Computational formulas 

Here we provide the explicit coordinate formulae for the conformal metric and the con- 
formal cubic form as introduced in Sect. 4.2. 
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The Lagrangian Grassmannian A is (locally) parametrized by 3 x 3 symmetric matrices 
(26), so that the cubic form C is given by 

C = det dU = dadcdf — da(dqY — (db^df + 2dbdpdq — dc{dpf'. 

Each tangent space to the Lagrangian Grassmannian (which can also be identified with 
3x3 symmetric matrices) carries a cubic cone C = consisting of rank two symmetric 
matrices. Its singular locus coincides with rank one matrices, and can be parametrized in 
the form 

1 /i A \ 

/i^ A/i (29) 
A A/i A2 / 

where A and /i are independent parameters. The restriction Cm^ of the cubic form C to 
the hypersurface M^, defined by the equation / = /(a, 6, c, p, g), is obtained by setting 

df = fada + fbdb + f^dc + fpdp + fgdq. (30) 

This results in 

Cm^ = fa{da)'^dc + fbdadbdc + fpdadcdp + fcda{dcY + fgdadcdq — da^dq)"^ — 

(31) 

fada{dbf - h{dbf - fp{dbfdp - fddbfdc - fq{dbfdq + 2dbdpdq - dc{dpf, 

here a,b,c,p,q are local coordinates on M^. The intersection of the singular locus (29) 
with the tangent space to is a rational normal curve with the affine parametrization 

7 = (1, Ai, A, /iA); 

notice that A and fi are no longer independent and satisfy, by virtue of (30), the quadratic 
relation 

= /a + /fe/i + /c/i^ + /pA + /gA/i, 

which coincides with the dispersion relation (8). The tangent variety of 7 is defined as 
7 + ^7' where prime denotes differentiation with respect to fi, and A is viewed as a function 
of fi specified by the dispersion relation. Explicitly, the tangent variety is given by 

(1, fi + t, j/ + 2tfx, X + t\', fi\ + t{\ + 

here A' is obtained by implicitly differentiating the dispersion relation with respect to fi, 

A' = ih + 2/,/i + /,A)/(2A - V - fp)- 

One can verify that, up to a conformal factor, there exists a unique quadratic form Qm^ 
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vanishing on the tangent variety of 7, 

= (4/2 + UQda^ + {8fJb + hfl + 2fJMdadb + {ifj, + f^)dadp 

Hfb + 4/a/c + fcf^ + faP,)dadc + + flf,)dadq 

H^fi + ^fafc + 2hfMdb^ + + 4/,/, + 2flf,)dhdp 

(32) 

+ (8/Jc + + /6/2)rfMc + (4/Jp + 2/,/, + 2Ufl)dhdq 

-(4/, + + + /p/,2)rfprfc - (4/, + 2fjg)dpdq 

+ (4/.' + /c/,')c?c2 + (4/J, + /3)rfcrfg - (4/, + 

(one needs to use the dispersion relation in order to verify the vanishing of Qm^ on the 
tangent variety). Representing Qm^ as a 5 x 5 symmetric matrix and calculating its 
determinant one obtains 

det = S{f^ + - fj^ - Ul - AfJ^Y = 3 A'; 

we point out that A = if and only if the corresponding dispersion relation (8) defines 
a reducible conic. Thus, under the assumption of the irreducibility, the conformal metric 
Qm^ is non-degenerate. The cubic form (31) and the conformal metric (32) satisfy the 
relations 

CijkQ^-' = 0, 

CjkrQ Clns ~\~ CljrQ Ckns ~\~ CklrQ Cjns 27A2 (^QjkQln ~l~ QljQkn ~l~ Q klQ jn) • 

Thus, the normalization C 3^/3 AC results in the identities (25). 



4.5 General aspects of irreducible G/(2)-structures 

Let be a projective space with coordinates ( A rational normal 

curve 7 of degree four can be parametrized in the form 7 = (1 : t : : : t^). Its bisecant 
variety is a cubic hypersurface defined by the equation 

n 1 2 

00 3j 

det I 

2 4 

OC OC 00 

notice that 7 can be recovered as a singular locus of its bisecant variety. The tangent 
variety of 7 is contained in a unique quadric hypersurface 



X 



V-4x^x=^ + 3(^2)2 = 0. 



The 'curved' analog of this picture is the following. Let be a 5-dimensional mani- 
fold with a field of rational normal curves of degree four specified in the projectivization 
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of each tangent space. One can choose a frame of 1-forms u^, u"^, u^, u'^, lJ' such that the 
equation of the bisecant variety of 7 takes the canonical form C = where 



C = det 



ijj 



such frame is defined up to a natural G'L(2)-equivalence, leading to the following structure 
equations: 



du^ = Au'^ + Auj^ + 'ip Ato^ + 3r] Au^ + 
du'-^ = n A + 2tp A uj^ + 2r] A LJ^ + 
du^ = n Au^ - (p Au'^ + 3ij Alj^ + r] Au^ + 
du^ = n A uj^ - 2(p A uj^ + iij A + .... 



(33) 



Here dots denote linear combinations of the terms A uo^ which are responsible for the 
'torsion'. The structure equations for the secondary forms Q,(p,ip,ri, restricted to the 
fiber u'' = 0, take the standard G'L(2)-form: 

dfl = 0, dr] = (p A T], dip = 21] A ip, dip = —ip A ip. 

Introducing uj = {u^ , u"^ , , , u^Y , one can represent (33) in a compact form 

du = {QE + ipEi + ipH + r]E2) Auj + T{uj, u), (34) 

where T represents the torsion, E is the identity matrix, and Ei,H,E2 define the 5- 
dimensional irreducible representation of SL{2): 
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We point out that bisecant surfaces and trisecant three-folds make perfect sense for ab- 
stract GL{2) structures. Recall that if the GL{2) structure comes from a hypersurface of 
the Lagrangian Grassmanian, it automatically possesses infinitely many bisecant surfaces 
parametrized by two arbitrary functions of one variable. For abstract GL{2) structures, 
the existence of bisecant surfaces imposes strong restrictions on the torsion T, which can 
be represented in a simple geometric form: 

Theorem 4 An abstract GL{2) structure possesses infinitely many bisecant surfaces if 
and only if, for any two vectors X and Y which belong to the rational normal curve, one 
has G{X,Y,T{X,Y)) = 0. 
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Proof: 

Let us look for bisecant surfaces in the form u = a a dR^ + (3 h dF? where i?^, are 
the characteristic coordinates, a and b are (column) vectors which belong to the rational 
normal curve 7, that is, 

a = (1 : a : : : a^)\ h = {I : h : h"" : : b^, 

and a, b, a, (3 are certain functions of R^,R^. Substituting this ansatz into (34) and 
using the identities 

Eia = a', ^ib = b', 

Ha = 2a- aa', iJb = 2b - 6b', 

E'aa = 4aa - a^a', Esb = 46b - 6^b', 

where 

a' = (0 : 1 : 2a : 3a2 : Aa^, b = (0 : 1 : 26 : 36^ : 46^)*, 

one arrives at 

a{da -aQ- 2ai^ - Aaat]) A dR^ + h{df3 - (3Q - 2f3ip - 4b(3r]) A dR'^+ 

aL'{ada - atp + aa^p + aa^r]) A rfi?^ + h\l3db - + f3bip + (3b'^r]) A dR^ = (35) 

2ap T(a, b) dR^ A dR^. 

Thus, for the consistency of these equations one has to require 

T(a, b) G span {a, b, a', b'}, 

which is equivalent to the statement of the theorem. Once this condition is satisfied, one 
obtains a system of four first order PDEs for a, b, a, (3 by collecting similar terms in (35). 
Up to reparametrizations i?* — *■ ip^{R^), the general solution of this system depends on 
two arbitrary functions of one variable. 

Theorem 4 establishes necessary conditions for the realizability of a GL{2) structure 
on a hypersurface of the Lagrangian Grassmannian. Further development of the gen- 
eral theory of abstract GL{2) structures is beyond the scope of this paper, and will be 
addressed elsewhere. 



5 Integrability in more than three dimensions 

The integrability of three-dimensional equations of the form (1) was defined as the exis- 
tence of n-component hydrodynamic reductions (2) parametrized by n arbitrary functions 
of a single variable. As demonstrated in [23, 24], this approach readily generalizes to any 
dimension: a d-dimensional PDE 

F(m,w)=0 (36) 
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for a function u of d independent variables x^, x is said to be integrable if it possesses 
infinitely many n-component hydrodynamic reductions parametrized by {d — 2)n arbitrary 
functions of a single variable. In this case equations (2) are replaced hj d—1 commuting 
(1 + l)-dimensional systems of hydrodynamic type. Among the known four-dimensional 
integrable examples one should primarily mention the 'first heavenly equation', 

UxyUzt UxtUzy 1; 

as well as its equivalent forms, 

"^tx ~l~ ^zy ~l~ Uxx'^yy ^xy ^ 

and 

Utt = UxyUzt — UxtUzy 1 

known as the 'second heavenly' and the 'Grant equation', respectively [47, 28]. It was 
demonstrated in [23, 24] that n-component reductions of these equations are parametrized 
by 2n arbitrary functions of a single variable. 

An interesting six- dimensional integrable generalization of the heavenly equation, 

Utt + Uzz + UtxUzy - UtyUzx = 0, 

arises in the context of sdif f{Y?) self-dual Yang- Mills equations [48]. Its n-component 
reductions are parametrized by 4n arbitrary functions of a single variable [24]. Notice 
that all these examples belong to the class of symplectic Monge-Ampere equations as 
introduced in Sect. 3. 

Although the general problem of classification of multi-dimensional integrable equa- 
tions can be approached in a similar way, the method of Sect. 2 leads to quite complicated 
analysis. One way to bypass lengthy calculations is based on the following simple idea: 
suppose we want to classify four-dimensional integrable equations of the form (36) for a 
function z, t). Let us look for travelling wave solutions in the form 

u{X, Y, Z) = u{x + at, y + /St, z + ^t). 

The substitution of this ansatz into (36) leads to a three-dimensional equation which must 
be integrable for any values of constants q;,/3,7. Since, in three dimensions, the integra- 
bility conditions are explicitly known, this provides strong restrictions on the original 
function F, which are necessary for the integrability. The philosophy of this approach is 
well familiar from the soliton theory: symmetry reductions of integrable systems must be 
themselves integrable. 

Thus, for the first heavenly equation, traveling wave solutions are governed by 

a{uxYUxz - uxxUyz) + liuxyUzz - uxzUyz) = 1, 

which is a three-dimensional symplectic Monge-Ampere equation. One can show that it 
is indeed integrable (in fact, linearizable) for any values of constants. 

We hope to return to the multi-dimensional case in subsequent publications. 
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6 S'2?(6)-action on the moduli space of integrable equa- 
tions 



In this section we investigate the action of the equivalence group Sp{6) on the moduh 
space of integrable equations of the dispersionless Hirota type (recall that both Sp{6) and 
the moduli space have coinciding dimensions equal to 21). Our main result states that 
this action has an open orbit. 

For a 3 X 3 symmetric matrix U = [uij], the Lie algebra of the group of conformal 
automorphisms of the cubic form det dU is spanned by 21 vector fields which generate 
the Lie algebra of the symplectic Lie group 5^(6): 
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Given a PDE of the form F{uii, uu, ^13, ^22, ^23, M33) = 0, we will look for its infinitesimal 
symmetries by solving the determining equation LxF\p=Q = 0, where X is a linear com- 
bination of the 21 vector fields presented above (notice that the answer may not coincide 
with the full algebra of Lie-point symmetries: we consider only those symmetries which 
belong to the equivalence group 5*^(6)). Below we list particular examples of integrable 
equations which possess symmetry algebras of different dimensions (it is worth noting that 
any two equations with different symmetry algebras are automatically non-equivalent). 
Example 1. The 2-dimensional linear wave equation, uu + U22 — U33 = 0, possesses nine 
infinitesimal symmetries: 

X12, X13, X23, X11 + X33, X22 + X33, J1 + J2 + J3, L12 — L21, L13 + L31, L23 + L32. 

One can prove that the existence of nine infinitesimal symmetries is necessary and suffi- 
cient for the linearizability of a general equation of the form (1), [43]. 
Example 2. The dKP equation, M22 — M13 + ^ufi = 0, possesses seven infinitesimal 
symmetries: 

X12, X23, X33, X13 + X22, X11 + L31, Ji -|- 2J2 + 3J3, 2L32 + L21. 

It is likely that there exist no integrable equations with eight symmetries. 

Example 3. The Boyer-Finley equation, un + U22 — e"^^ = 0, possesses six infinitesimal 

symmetries: 

X12, X13, X23, Xii — X22, J1 + J2 + 2X33, L12 — L21. 

Example 4. The degeneration of the dispersionless Hirota equation, U12 — M13 — e"^^ = 0, 
possesses five infinitesimal symmetries: 

Xii, X22, X33, J1+X23, X12 + X13. 

Example 5. The equation e"" -|- e"^^ — e"^^ = possesses four infinitesimal symmetries: 

X12, X13, X23, Xii + X22 + X33. 
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Example 6. The dispersionless Hirota-type equation for the BKP hierarchy, e"" + e"^^ + 
g«23 — gMi3+wi3+ii23^ possesses three infinitesimal symmetries: 

Xii, X22, X33. 

Remark. We point out that the existence of 'many' symmetries is not related to the 
integrability: for instance, both equations (20), which are not integrable, possess 8- 
dimensional symmetry algebras isomorphic to SL{3,R) and SU{3,R), respectively [4]. 
Thus, the equation Hess u = 1 possesses eight infinitesimal symmetries: 

L12, Li3; L21, L23, L31, L32, Ji — J2, Ji — J3. 

The equation Hess u = un + ^22 + u^s also possesses eight infinitesimal symmetries: 

X22— Xii, X33 — Xii, P1 + X23 P2+X13, P3 + X12, L12 + L21, L13 + L31, L23 + L32. 

The main result of this section is the following 

Theorem 5 The action of the equivalence group Sp{6) on the moduli space of integrable 
equations of the dispersionless Hirota type has an open orbit. 

This fact is, in a sense, surprising: it establishes the existence of a 'universal' equa- 
tion with no symmetries, which generates an open part of the moduli space under the 
action of 5*^(6). In particular, one should be able to obtain all equations with non-trivial 
symmetries by taking appropriate degenerations of this universal equation. 

Proof of Theorem 2 

The main idea of the proof is to prolong the 21 infinitesimal generators Xn — P3 to 
the moduli space of solutions of the involutive system (16). We point out that, since third 
order derivatives of / are explicitly known, this moduli space can be identified with the 
values of / and its partial derivatives /j, fij up to second order (21 parameters altogether). 
The prolongation can be calculated as follows: 

(1) Following the standard notation adopted in the symmetry analysis of differential 
equations [31, 44], we introduce the variables 

X^=Uu, X^=Ui2, X^ = Ui3, X^=U22, = U23, U = U33, 

and rewrite the above 21 generators in the form 

here C,^ and rj are certain functions of x and u. In this notation, a dispersionless Hirota- 
type equation is represented in the form u = u{x^, x^) (the function u is denoted by / 
in Sect. 2). 

(2) Prolong infinitesimal generators to the second order jet space with coordinates u, Ui, Uij, 

d d ^ d ^ d 

C—+V— + C^— + C^3 



dx^ du dm dm 
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where d and Qj are calculated according to the standard prolongation formulae 

= A(r/) - UkD.ie), = D,Q - u,kD,{ey, 

here Di are the operators of total differentiation. 

(3) To eliminate the ^-terms, subtract the linear combination of total derivatives ^^Di 
from the prolonged operators. It is sufficient to keep only the following terms in Di: 

d d d d 

Di 7^ J ~\~ Ui~z h '^ijT^ I" ^iiTcTI i 
OX* OU OUj OUjk 

notice that, since Uijk are explicit functions of lower order derivatives, the resulting opera- 
tors will be well-defined on the 21-dimensional space with coordinates u, Ui, Uij. Although 
these operators will depend on the variables as parameters (indeed, the isomorphism of 
the moduli space with the space u, Ui, Uij depends on the choice of a point in the x-space), 
all algebraic properties of these operators will be x-independent. 

(4) Finally, the dimension of the maximal S'p(6)-orbit equals the rank of the 21 x 21 
matrix of coefficients of these operators. It remains to point out that this rank equals 21 
for any 'random' choice of numerical values for x^,u,Ui,Uij (however, it equals 21 — r for 
any example with r symmetries). 

7 Concluding remarks 

We would like to formulate a list of natural questions which are left beyond the scope of 
this paper. 

— It was demonstrated that the integrability of equations of the form (1) can be inter- 
preted as the existence of 'sufficiently many' trisecant sub manifolds of the corresponding 
hypersurfaces in the Lagrangian Grassmannian. It would be desirable to have a tensor 
characterization of this condition in terms of the metric Qij and the cubic form Cijk in- 
troduced in Sect. 4.2. Presumably, these conditions will involve the Weyl tensor of the 
conformal metric Qij (which is not conformally flat in general). Alternatively, one may try 
to express the integrability conditions using the structure equations of the corresponding 
GL(2)-structure. 

— We have shown that any hypersurface of the Lagrangian Grassmannian possesses an in- 
trinsic G'L(2)-structure. In would be desirable to develop a general theory of hypersurfaces 
of the Lagrangian Grassmannian, as well as the theory of abstract GL(2)-structures, and 
to establish the embedding theorem. Notice that the theory of curves in the Lagrangian 
Grassmannian received some attention, see [45, 42, 55]. 

— One of our main observations is the existence of a universal 'master-equation' which 
generates all other integrable equations of the dispersionless Hirota type via various (sin- 
gular) limits. It seems to be very important to understand the reduction procedure 
geometrically, and to obtain the master-equation explicitly in terms of the appropriate 
special functions (one shouldn't expect simple formulas since even particular examples 
contain transcendental functions such as the elliptic theta-functions and solutions to the 
Chazy equation). 
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To some extent, the picture resembles the situation with the classification of (non- 
holonomic) eyelids of Dupin in E"^: as demonstrated in [46], up to the action of the Lie 
sphere group SO {5, 2), there exists a unique representative, namely, the stereographic 
projection of the isoparametric hypersurface of Cartan. Coinciding dimensions of the 
equivalence groups, dim Sp{6)=dim 5*0(5, 2) =21, add to this similarity. Notice, however, 
that any nonholonomic Dupin hypersurface in E"^ is Lie-homogeneous, while our generic 
example possesses no continuous S'p(6)-symmetries. 
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